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posed by Steven R. White (Phys. Rev. Lett. 69 (1993) 3844). We find a massless 
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T^-j- ■ point at 5c = 0.25 ± 0.01. We also find the edge states in the region 5 < 5c under 
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the open boundary condition, which disappear in the region 5 > 5c. At the massless 
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C ' point, the spin wave velocity Vg is 3.66 it 0.10 and the central charge c is 1.0 it 0.15. 
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The antiferromagnetic 5* = 1 Heisenberg chain has been the subject of a large amount 
of interest in the last decade since Haldane argued the difference of the low energy physics 
between integer and half-integer spin ^^. The key concept of the Haldane conjecture is 
described in terms of the 0(3) nonlinear sigma model with the topological term and he 
concluded that the topological term is ineffective and the spectrum is massive for integer 
spin, while the topological term suppresses the mass-generation and excitation spectrum is 
massless for half-integer spin. For the S = 1/2 Heisenberg model, it is established that the 
excitation spectrum is gapless by the Bethe-Ansatz solution. The existence of excitation gap 
in the case of 5 = 1 was verified from the theoretical |P,^ , numerical |^|-@ , and experimental 



results 10-13 



Some years later, Haldane and Affleck |[l4|-p!8[] considered the antiferromagnetic (AF) 
spin chain with bond alternation S, 



Ns-l 



n= Yl [i-i-iy5)s,-s,+,. (1) 



i=l 



When this Hamiltonian is mapped to the continuum model, i.e., the 0(3) nonlinear sigma 
model, the topological angle 6, which is the coefficient of the Berry phase term in the action, 
is given by ^ = 27iS{l — 6). According to the argument by Haldane and Affieck, only 
in the case where 6 = it x (odd integer), is the system in the critical state with massless 
excitations. Hence as the bond alternation 6 changes from —1 to 1, the system experiences 
critical (massless) point 25" times. Applied to the S = 1/2 case, this general discussion 
is consistent with Cross-Fisher's argument of the spin-Peierls transition [|1^. However the 
situation is nontrivial in the case of S* = 1. Without the bond alternation 6, the system is 
gapful, i.e., possesses the Haldane gap, and the above discussion predicts that the massless 
point appears at some 6 = 6c ^ 0. However this argument contains several approximations 
and assumptions, such as taking the large S limit, and it is important to study the model 
(|I|) by another approach and to give an independent check of the above scenario. 

Varying-(5 problem in the model (1) can be attacked from the viewpoint of the topological 
order proposed by den-Nijs and Rommelse ||20| and Tasaki |2l|. This "hidden" order is 



measured by the string correlation function 

9{^,J) = {SnU e^'^^ns;). (2) 

\k=i+i / 

It is widely believed that this correlation function is long-ranged at the Heisenberg point 



(5 = from numerical evidences p2| , |23[] . On the other hand, an elementary calculation 
shows that the string correlation function is exactly zero at the completely dimerized points 
6 = ±1. Hence we expect that phase transitions corresponding to the hidden symmetry 
breaking should occur twice; once in the interval between 6 = —1 and 6 = and once 
more between 6 = and 6 = 1. However, these arguments cannot determine whether 
the phase transitions are discontinuous or continuous. Thus the existence of the massless 
points as predicted by Haldane and Affleck is still an open problem and direct numerical 
study is needed. In this paper we investigate the model ([^) with S = 1 numerically for 
the region < 6 < 1, using exact diagonalization and the density matrix renormalization 
group (DMRG) method. This is sufflcient since reversing the sign of 6 merely corresponds to 
interchanging the roles of the two sublattices, and hence has no effects on the bulk properties. 
It does however have an intrinsic boundary effect for open chains on which we will remark 
later. 

DMRG, which was recently proposed by White, is a new algorithm for real-space renor- 
malization group [^3l -|26[|. The key point lays in the method of the truncation procedure of 
the Hilbert space. In conventional real space renormalization group calculation, one keeps 
the lowest-lying eigenstates of the block Hamiltonian in forming a new effective Hamilto- 
nian. In the DMRG method, on the other hand, one keeps the most important eigenstates 
of the block density matrix, obtained from diagonalizing the Hamiltonian of a larger section 
of the lattice which includes the block. The above change brings high accuracy and the new 
algorithm is already applied to the S=l antiferromagnetic Heisenberg spin chain successfully 



6| . For details of the DMRG algorithm, readers are referred to ref . [^ 
Using the DMRG method, we calculated the excitation gap of systems of large length 
up to about two hundred sites and extrapolated them to the thermodynamic limit. The 



numbers of state kept per block in our calculations are 27, 36, 48, and 81. As a benchmark 
test we calculated the ground state energy density oi S = 1/2 Heisenberg model. Comparing 
with the exact Bethe-Ansatz solution, we find that the relative error is about 5 x 10~^. We 
also compared those of DMRG with that of exact diagonalization in the case of 5* = 1 and 
site number A^^ = 12. The two results were found to agree well within a relative error of 
7 X 10^^ 

Figure 1 shows the excitation spectrum of S" = 1 spin chain versus dimer strength. The 
data of diagonalization (points plotted with triangles, squares, empty and filled circles and 
empty and filled diamonds) shows the gap between the singlet ground state and the first 
excited (triplet) state under periodic boundary condition. The results of diagonalization 
show that the gap has a minimum around 6 ~ 0.25, which has been already reported in ref. 



27| . The minimum value of the gap becomes smaller with increasing system size. In the 
large A^^ limit, there are two possibilities i.e. the minimum value of the gap tends to zero or 
saturates at a finite value. Limited system sizes in the diagonalization, however, disenables 
us to conclude which is realized in the thermodynamic limit. 

The data of DMRG calculation (represented by the double circles and crosses), which 
are extrapolated values {Ng -^ +oo), show the gap between the ground state and the first 
excited state, and that between the ground state and the second excited state under the 
open boundary condition. The gap behaviors are different between the region 6 < 0.25 and 
6 > 0.25. In the region 6 < 0.25, the first excited state is triplet {Stot = 1) and the second 
one is quintuplet (S'tot = 2) while both states are triplet in the region 6 > 0.25. The first 
excited state in the region 6 < 0.25, approaches the singlet ground state exponentially with 
increasing system size. These triplet excitations have been investigated at the Heisenberg 
point {6 = 0) and are known to be the edge excitations near the open boundary of the spin 
chain [^,^]. For the moment, we focus upon the gap between ground state and the second 
excited state in the region 6 < 0.25, because the above "edge excitations" seem to have little 
effects upon the bulk properties of the spin chain. Later we will discuss the edge states, 
especially their relevance to the hidden symmetry breaking. 



The size dependence of the gap between the ground state and the second excited state in 
the region 6 < 0.25 is fitted well with 1/N^- The gap between the ground state and the first 
and second excited states in the region 6 > 0.25 is also fitted well with l/N^. At 6 = 0.25, 
size dependence of the gap behaves as l/Ng. In contrast to the exact diagonalization method, 
DMRG gives a clear evidence for the existence of a massless point at 6c = 0.25±0.01. 

We now turn to the critical properties of the massless point. It is plausible to assume 
that at the massless point, the system has conformal invariance with respect to the low 
energy properties. Actually, we observe that the size dependence of the excitation gap is 
proportional to 1/A^s, which is consistent with the finite size scaling of the conformal field 
theory (CFT). CFT gives the expression for the size dependence of ground state energy 

under the open boundary condition. Here Eg,s{Ns), e, /, and c are the ground state energy of 
finite size system, ground state energy density of the infinite system, surface energy and the 
central charge. We estimate Vg = 3.66 ± 0.10 at 6 = 6c from the excitation gap. Combining 
the value of Vg with the results for the ground state energy, we estimate c = 1.0 ± 0.15. 

Affleck and Haldane predicted 6c = 0.5, which deviates from our result 6c = 0.25 ± 0.01. 
We can attribute this discrepancy to the mapping of the quantum spin chain onto the 0(3) 
sigma model, which is justified in the large S limit. However their prediction is correct with 
regards to the existence of a massless point. 

The existence of edge states or four-fold degeneracy of the ground state in the open 
chain was discussed theoretically 0,|,^,^] and was related to the long range topological 



order and hidden Z2 x Z2 symmetry breaking. So we conclude that at the massless point a 
continuous phase transition with the hidden Z2 x Z2 symmetry breaking occurs. Our result 



(c ~ 1) is in accord with the general belief concerning the Z2 x Z2 symmetry breaking ||3^ 
that the massless phase boundary which corresponds to a full Z2 x Z2 symmetry breaking 
belongs to the same universality class as the Gaussian model. 



We note in passing the consequence of our choice of the sign of S; had we taken negative S, 
we would have found that the ground state degeneracy is nine- fold in the region —1<S<6~ 
(due to the formation of two localized S* = 1 objects at the open ends), while it is still four- 
fold in the region S~ < 6 < 0, where 6~ = —6c ~ —0.25. 

In summary, we presented numerical results on the S = 1 antiferromagnetic spin chain 
with bond alternation and concluded that the continuous phase transition occurs at 6c = 
0.25 ± 0.01 with the hidden Z2 x Z2 symmetry breaking. The higher S cases are left for 
future investigations. 
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FIGURES 
FIG. 1. Excitation gap between the ground state and excited states versus the dimer strength 

6. We show here two kinds of numerical results. Points plotted with triangles, squares, circles 

(empty and filled) and diamonds (empty and filled) are the gap between the ground state and the 

first excited state obtained by the exact diagonalization under the periodic boundary condition for 

system sizes Ng = 4, 6, 8, 10, 12, and 14. As the system becomes larger, the minimum of the gap 

shows up around S ~ 0.25. Double circles/crosses represent the gap between the ground state and 

the first/second excited state obtained by the DMRG method (the infinite lattice method) under 

the open boundary condition. We find a massless point around 6c = 0.25 it 0.01, and also edge 

states (Kennedy triplet) in the region S < 5c- 



